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ABSTRACT 


One of the most serious problems of modern navigation is the accident rate that 
occurs due to inept or belated maneuvering of ships. As a result of accidents in the world, 
more than 200 ships sink every year and every fourth receives significant damage. 

Full-scale tests show that the stopping distance of large-tonnage ships turn out to be 
much less permissible, and shipbuilders are able to significantly reduce the astern power of 
such ships, making them cheaper at the expense of safety. 

The low accuracy of inertial-braking characteristics is mainly due to unqualified field 
tests. Analysis of graphs and tables based on the results of such tests show that the spread 
in the values of inertial-braking characteristics for ships of the same type reaches 30%, and 
in some cases even more. In many tables and graphs, inertial-braking characteristics are 
expressed in relative values and are not suitable for direct use when maneuvering a ship. 
Finally, even when graphical and/or tabular maneuvering information is available on the 
navigating bridge, it is difficult to use it when maneuvering a ship at night. 

The research carried out by the author results in: 

- creation of an alternative computational method for determining the iner- 

tial-braking characteristics of the ship, suitable for use on any on-board computer; 

- development of an improved methodology for calculating the path and time of 

acceleration and braking of the ship in various ahead motion modes; 

- development of a methodology for taking into account the influence of a passing 

and opponent current on the length of the stopping distance of the ship; 

- development of methods for solving applied problems, ensuring a decrease in the 

accident rate of ships during maneuvering. 

The obtained methods include the development of theoretical foundations, mathe- 
matical models and comparison of the calculated inertial-braking characteristics of ships 
with the data of a full-scale experiment. For the first time, to derive the calculated formulas 
for the time and stopping distance, theorems are used on the change in the momentum 
and kinetic energy during accelerated and decelerated motion of the ship. In the course of 
the study, the problems of calculating and formalizing the inertial-braking characteristics 
of the ship are being comprehensively solved. For the first time, the hypothesis that the 
nature of the change in the thrust force of the propeller during reverse can be approximated 
by linear equations has been substantiated and confirmed. The general results are used to 
calculate the inertial-braking characteristics of specific ships. 
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tics, full-scale tests, propeller thrust, resistance to ship movement. 
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INTRODUCTION 


One of the most serious problems of modern navigation is the accident rate 
that occurs due to inept or belated maneuvering of ships.As a result of accidents 
in the world, more than 200 ships sink every year and every fourth receives signi- 
ficant damage. Numerous cases of collisions of ships, bulkheads on each other and 
on shore facilities, groundings, etc. also indicate that a significant defect of some 
ship complexes is the long duration of the reversing process and an unacceptably 
long stopping distance. 

On December 5, 2002, the Maritime Safety Committee of the International 
Maritime Organization adopted standards for the maneuverability of trans- 
port ships, which would increase safety at sea and protect the marine environ- 
ment. However, according to these criteria, the braking distance when braking 
from full ahead to full astern is allowed up to 15-20 hull lengths, depending 
on the ship’s power-to-weight per tonne of displacement. Full-scale tests of an 
oil tanker with a displacement of 156,230 tons and a tanker with a displace- 
ment of 125,273 tons showed stopping distances for these ships, respectively, 
10 and 7 hull lengths. Consequently, shipbuilders were able to significantly re- 
duce the reverse power of such ships, making them cheaper at the expense of 
reduced safety. 

On the navigating bridges of many ships, information about their iner- 
tial-braking characteristics is either completely absent, or it is not presented fully 
and accurately enough. The low accuracy of inertial-braking characteristics is 
mainly due to unqualified field tests. Analysis of graphs and tables based on the 
results of such tests has been shown that the spread in the values of inertial-brak- 
ing characteristics for ships of the same type reaches 30 %, and in some cases even 
more. In many tables and graphs, inertial-braking characteristics are expressed in 
relative terms, for example, in hull lengths, so they are generally not suitable for 
direct use when maneuvering a ship. 

Finally, even when graphical and/or tabular maneuvering information is 
available on the navigating bridge, it is difficult to use it when maneuvering a ship 
at night. Such information could be stored in an onboard computer connected 
to lag sensors, gyrocompass and other navigation devices, and continuously dis- 
played on the monitor screen. 

The research carried out by the author resulted in: 

- creation of an alternative computational method for determining the iner- 

tial-braking characteristics of the ship, suitable for use on any on-board 
computer; 
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- development of an improved methodology for calculating the path and 
time of acceleration and braking of the ship in various ahead motion 
modes; 

- development of a methodology for taking into account the influence of 
a passing and opponent current on the length of the stopping distance 
of the ship; 

- development of methods for solving applied problems, ensuring a de- 
crease in the accident rate of ships during maneuvering. 

The research methodology includes the development of theoretical founda- 
tions, mathematical models and comparison of the calculated inertial-braking 
characteristics of ships with the data of a full-scale experiment. In the theoretical 
conclusions, the theorems of changing the momentum and kinetic energy have 
been used. The hydrodynamic drag force has been taken to be proportional to 
the ship’s speed to any positive degree, and the thrust force of the propeller has 
been approximated by linear equations. The work of the forces of water resistance 
and the thrust of the propeller has been equated, respectively, with the work of 
their mean integral and arithmetic mean values. When developing the theoretical 
foundations, the results of numerous domestic and foreign studies were analyzed 
and summarized. 

In this work, the tasks of calculating and formalizing the inertial-braking 
characteristics of the ship are comprehensively solved. For the first time, the hy- 
pothesis that the nature of the change in the thrust force of the propeller during 
reverse can be approximated by linear equations has been substantiated and 
confirmed. This hypothesis serves as the basis for the derivation of the calculation 
formulas. For the first time, to derive the calculated formulas for the time and 
stopping distance, theorems has been used on the change in the momentum and 
kinetic energy during accelerated and decelerated motion of the ship. Calculation 
formulas allow to continuously obtain the necessary information when changing 
the propeller speed and the speed of the ship. The general results obtained in this 
work are applied to calculate the inertial-braking characteristics of specific ships. 
For the first time, a methodology has been developed to take into account the 
influence of a tail and opponent current on the stopping distance of a ship. 





Chapter 1 
REVIEW AND ANALYSIS 
OF LITERATURE 


In this chapter, the analysis and structuring of the works of recognized domestic 
and foreign scientists, theorists and practitioners involved in the development of 
methods for calculating the inertial-braking characteristics of ships and working 
in various directions is carried out. The author also points out a number of 
serious inaccuracies made in the educational literature, which lead to incorrect 
perception and interpretation by inexperienced navigators of the ship behavior 
when maneuvering, and, consequently, to the emergence of considerable errors 
in the calculations of the maneuvering characteristics of the ship. 

Keywords: propeller, reverse, stopping distance, maneuverable elements, method 
of calculating maneuverable characteristics. 


The current period of development of maritime transport is characterized by 
the improvement of transportation technology and specialization of ships by types 
of cargo and directions of transportation [1, 2]. The specialization led to a sharp 
increase in the size and speed of ships, and also significantly reduced the parking 
time under cargo operations, increasing the time spent at sea. In such conditions, 
the requirements for the design of the maneuverable characteristics of ships are 
radically changing. 

Unfortunately, the possibilities of the designer when designing the maneu- 
verable characteristics of ships are very limited. The main dimensions and shapes 
of the hull lines are selected in advance for reasons of propulsion, seaworthiness, 
stability, etc. Only the propulsion unit and the steering device remain at its dis- 
posal. The consequence of this is the long duration of the reversal process in some 
high-speed ships [3, 4]. 

Reverse is a maneuver to change the direction of movement of the ship by 
changing the direction of action of the thrust force of the propeller. It is one of 
the most difficult and therefore the most critical operating modes of the ship 
complex [5]. 

Intuitively, it seems that with an increase in engine power astern, supplied 
to the propeller, the distance and braking time of the ship will be reduced when 
reversing from ahead to astern. However, this dependence does not exist in all 
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stages of the reverse. A stream of water with complex particle movement is created 
behind the stern of a forward ship. This movement can be simply represented 
as a spiral [5, 6]. When the fixed pitch propeller (FPP) is reversed from ahead 
to astern, energy has to be expended on turning the flow, both in the axial and 
circumferential directions. In this case, the energy consumption turns out to be 
much higher than those calculated theoretically. In addition, if the reversal is per- 
formed too quickly, then the flow reversal process is accompanied by cavitation, 
which contributes to additional energy dissipation, hence, an increase in the time 
and length of the stopping distance before the ship stops [7]. 

In order for the engine power to be used with maximum efficiency, it is ne- 
cessary to select the optimal flow reversal mode, which would exclude eddies and 
cavitation. Due to the lack of reliable theoretical methods, it is not yet possible to 
calculate this regime [8]. 

The design of controllable pitch propellers (CPP) allows, during reverse, 
to change only the axial direction of flow without changing the circumferential, 
which means that the change in the direction of the thrust force can be carried out 
without cavitation and with a greater effect than the FPP reverse. However, due 
to the complexity of the design and some operational inconveniences (continuous 
rotation of the propeller even in the neutral position), CPPs are not yet widely 
used [9, 10]. 

After reversing the flow in the opposite direction, further increases in power 
will shorten the distance and time to a complete stop of the ship. But after a certain 
limit, further increase in power becomes economically unprofitable, because the 
reduction in stopping distance will be disproportionately small compared to the 
energy consumption [10]. 

In view of the above difficulties, only approximate methods have been de- 
veloped for calculating the hydrodynamic forces arising from the reverse. Experi- 
ment remains the main source of information about these forces. 

Research of inertial-braking characteristics of a ship is carried out in two 
directions [11]. 

The first of the directions was developed in detail by V. Nebesnov [12] and 
A. Hoffman [13]. It is based on a thorough study of the individual stages of tran- 
sient processes and the choice of optimal operating modes for ship propulsion 
systems. To achieve this goal, the authors have developed a technique for the joint 
solution of differential equations of the balance of forces and moments of ship 
propulsion systems. 

The second direction of research is to refuse to take into account the inertia of 
the rotating parts of the complex due to their relative smallness. For example, the 
time for the propeller to accelerate to the specified RPM is negligible compared to 
the time it takes the ship to accelerate to the speed corresponding to these RPM. 
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V. Nebesnov wrote as follows: «...the time constant of the engine is ten times less 
than the time constant of a displacement ship, therefore it can be assumed that the 
angular speed of the engine shaft at the initial stage of its acceleration or braking 
almost instantly reaches a certain intermediate value, and then changes slowly, in 
accordance with the change ship speed» [5]. 

Most of the domestic and foreign authors who took part in the development 
of methods for calculating the inertial-braking characteristics of ships used the 
second direction of research. These include V. Bakaev and V. Lavrentiev [12], 
M. Grechin [14-16], M. Leskov, A. Oganov and S. Kurguzov [17], S. Demina [18], 
V. Pavlenko [11], N. Solarev [19]. Foreign authors include Tani [2, 20], Jaegor and 
Jourdain [21], Hewins, Ruz, Chase and others [22]. 

V. Bakaev and V. Lavrentiev proposed to rebuild the results of model tests of 
propellers carried out by Nordstrom in the Gothenburg experimental basin. They 
calculated the universal propeller thrust, torque and propeller gait factors. These 
coefficients remain finite at any value of the propeller speed. In their work [12], cal- 
culation schemes are given for determining the path and time of the ship’s braking 
under the action of the propeller. In these schemes, it is assumed that the reverse of 
the main engine is carried out instantly, and the calculation of forces is carried out 
by step-by-step numerical integration, which requires cumbersome calculations. 

M. Leskov, A. Oganov and S. Kurguzov developed a method for determining 
the inertia of a ship using a universal table. To enter the table, it is necessary to 
measure the speed drop time during full-scale braking. The inertial characteristics 
are estimated using this method in a relatively simple way and without large ex- 
penditures of operating time. A significant drawback of the technique is the mea- 
surement of the speed drop using the ship log, which, as is known, has significant 
inertia. The use of the absolute lag will not take into account the current speed in 
the area of full-scale braking. 

The works of the Japanese researcher Tani [2, 20] are devoted to the study of 
the inertial-braking characteristics of supertankers. Based on the results of full- 
scale and model tests, the author comes to the conclusion that the thrust force 
of the propeller during reverse is equal, according to the work done, to its value 
in the mooring mode, corrected by a constant coefficient equal to 0.925. This 
approach simplified the derivation of design formulas, but introduced significant 
errors in the determined values, especially for ships of average tonnage. 

V. Pavlenko [11] carried out a detailed analysis of the factors influencing the 
inertial-braking characteristics of river ships, developed algorithms for calculating 
these characteristics on a computer, and analyzed some modern means of emer- 
gency braking. 

A. Maltsev [23] developed a detailed classification of the maneuvering ele- 
ments of ships and the algorithms for maneuvering when they diverge. 
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The listed and other works related to the second direction of research are 
based on the assumption that the force of resistance to the movement of the ship 
varies in proportion to the square of the speed. The entire reverse process is di- 
vided into three periods: 

1) movement of the ship from the moment the command is given until the 

moment when the fuel (steam) supply to the main engine is stopped; 

2) movement of the ship from the moment of stopping the supply of fuel 

(steam) to the main engine until starting it to astern; 

3) movement of the ship from the moment the engine is started astern until 

it stops completely. 

In the first period (5-7 s.), it is assumed that the speed of the ship and the rota- 
tional speed of the propeller are kept constant and equal to their initial values during 
steady motion. In the second period, it is assumed that the braking of the ship oc- 
curs only due to the force of resistance to motion and that the average value of the 
thrust force of the propeller is equal to zero. In the third period, it is assumed that 
the thrust force of the propeller reaches a predetermined value and then remains 
constant until the ship comes to a complete stop. In this case, the time of the second 
period is either calculated analytically or determined from field observations. 

Such an approximation of the thrust force of the propeller is permissible only 
with a rough estimate of the reversing characteristics of a diesel ship, but for turbo 
ships and ships with a pitch control propeller, the assumption of a zero value of 
the thrust force in the second period of reverse can introduce significant errors. 

On turbine ships, the propeller braking during reverse begins simultaneously 
with the closing of steam to the forward turbine and the opening of steam to the 
reverse turbine. The change in the propeller thrust force can be approximated by 
two linear equations: from the initial speed of the ship to the speed of the be- 
ginning of the reverse and from the speed of the beginning of the reverse to the 
complete stop of the ship, where the mean value of the propeller thrust will be 
numerically equal to the thrust in the mooring mode [14-16]. 

Other sources of systematic errors are the assumptions about the quadratic 
nature of the change in water resistance from the ship's speed, which is valid only 
for small values of the Froude number, as well as the assumption of the value of the 
added water mass constant and equal to 10 % of the mass displacement of the ship. 

Therefore, for the design of the main power plants and propellers, taking into 
account the given inertial-braking characteristics, as well as for the production of 
verification calculations of the characteristics of already built ships, a sufficiently 
justified and universal technique is needed to limit the influence of the noted 
disadvantages [3, 4]. 

Calculation data of inertial-braking characteristics can serve as a basis for pre- 
paring «Information to the captain about the maneuverable elements of the ship». 
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Such «Information» should include detailed graphs of inertial-braking characte- 
ristics, acceleration and braking characteristics, tables and diagrams of turnability 
elements, possibly in electronic form; recommendations for reversing, for ma- 
neuvering in narrow areas, in shallow water, in ice, in a storm, in conditions of 
low visibility and under other difficult sailing conditions, taking into account the 
characteristics of a particular ship. It is difficult to overestimate the importance 
and usefulness of such «Information» for training navigators in competent ma- 
nagement of their ship. 

In the «Information to the captain about the maneuvering elements of the 
ship», at the request of IMO, information on the acceleration and braking of the 
ship must be included. Unfortunately, as a rule, there is no such information on 
ships. In addition, in the existing literature on this issue [7, 24, 25], calculation 
methods are given, according to which it is impossible to obtain the values of 
the time and distance traveled by ships during acceleration and braking with 
sufficient accuracy for practice. Meanwhile, these parameters are important 
for ensuring the safety of maneuvering. For this reason, the search for suitable 
mathematical models of the ship's motion during acceleration and braking are 
very relevant. 

The researches [24, 25] concentrate the above, as well as other disadvantages 
of the existing calculation methods for determining the inertial-braking charac- 
teristics of ships. So, in order to simplify the integration of differential equations 
during acceleration and braking of the ship, the propeller thrust in them is as- 
sumed to be constant. However, this assumption does not correspond to reality, 
since the propeller thrust will change in accordance with the law of change in its 
relative pitch, depending on the aheadspeed of the ship. 

In the research [24], in the differential equation for ship acceleration, the 
following designations are adopted: 


dV 
m,——=-kV* +P, ia 
* dt : vey 
where m, - mass of the ship, taking into account the added mass of water; k - coef- 
ficient of proportionality; P, - propeller thrust on the ahead course (designations 
given in the textbook). 

In this differential equation, in addition to the assumption that the propeller 
thrust is constant, i.e. P,=const, one more assumption is made that in steady 
motion the resistance force is equal to the propeller thrust force, i.e. when the 
acceleration becomes zero, then: 

kV =P 


e? 
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which corresponds to reality, but only for a steady motion. However, the authors 
of the research on this basis convert the differential equation to the form: 


WV _y 


—=k(V?-V’), 1.2 
me at ) os 


st 


where V,, - value of the steady speed of the ship; V - current speed value. 

Further, the authors integrate the thus transformed differential equation and 
obtain calculation formulas for determining the path length and acceleration time 
of the ship. 

The main delusion of the authors is that they combined two absolutely in- 
compatible laws: the law of steady motion with the law of accelerated motion. This 
error is clearly visible if in their new differential equation (1.2) represents V=0, 
when acceleration starts from a stationary state of the ship, then the differential 
equation takes the form: 


dV 
<= Ry?, 
mM, dt st 


so the hull resistance to the movement of the ship was not negative, but positive. 
The propeller thrust force disappeared altogether. However, it is known that when 
the ship accelerates from a given speed (or from zero speed) to a steady speed, 
the inertia force will be equal to the difference between the propeller thrust 
force and the hull resistance force. Therefore, the calculated values of the path 
length and the acceleration time of the ship according to the formulas of the 
textbook [24] will be an order of magnitude larger and completely unsuitable for 
practical calculations. 

Further, in paragraph 5.4 [24], a practical calculation of the inertial-braking 
characteristics is given. In particular, the calculation of the proportionality coef- 
ficient «k» in the resistance formula R= kV’. 

Based on the statistical analysis of the results of a limited number of field 
observations, the authors propose the following empirical formula: 


k= 5880-+0.6540,/2. (1.3) 


where Q - area of the wetted surface; B/d — the ratio of the ship’s breadth to the 
mean draft (designations adopted in the textbook). 

In the theory of the ship, there is an approved general formula for the re- 
sistance of water to the movement of a ship. In accordance with the theory of 
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similarity, the resistance of the water to the movement of the ship is calculated 
using this general formula: 


R=Emy?, from where beh, (1.4) 


where € - coefficient of total resistance of the ship; p — mass density of water. 

Paragraph 5.4 of the research [24] also provides an example of calculating the 
inertial-braking characteristics for a container ship and, in particular, calculating 
the coefficient «k» according to the formula (1.3): 


k =5880+ 0.654: 5008 = = 11337. 


For comparison, calculate the coefficient «k» by the formula (1.4): 


_ 0.0035 -1020-5008 
- 2 


k = 8939. 


The difference between these calculations is more than 20 %, which will in- 
troduce a significant error in the calculated values. 

To calculate the thrust of the propeller during reverse in the textbook [24], 
the empirical formula is again introduced: 


2 
P= Pag( 1-17} (1.5) 


i 


where V - current value of the speed during active braking; V; - speed at the 
moment of propeller reversal (initial speed of active braking); Pax — maximum 
thrust force of the propeller, which is reached at the moment the ship stops relative 
to the water. 

In other words, the thrust force of the propeller appears only at the moment 
of the beginning of the reverse, and up to this moment the propeller rotating 
in the hydroturbine mode did not create either thrust force or resistance to the 
ship movement. 

The thrust force of the propeller is measured fairly reliably with a thrust 
gauge. According to such measurements [20, 22, 26-29] in the initial period of 
reversal, after closing the fuel (steam) to the main engine, the propeller still creates 
a positive thrust due to the inertia of the rotating masses of the engine and pro- 
peller shaft. As the speed of movement decreases, the propeller thrust decreases 
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and reaches the maximum negative value at the moment the reverse begins. Then 
the gauge value fluctuates around the average value, approximately equal to its 
value in the mooring mode. 

Fig. 1.1 schematically shows the change in the propeller thrust force during 
reverse from ahead to astern. The bold curve characterizes the change in thrust ac- 
cording to the thrust gauge data, and the thin curve - according to formula (1.5). 
According to the calculated curve, the thrust force of the propeller during reverse 
occurs only at the moment of reverse and reaches its maximum value when the 
ship is stopped relative to the water. 

Proposed in formula (1.5), the law of change in the thrust during reverse has 
nothing to do with the actual thrust measured with the thrust meter. 


P 





Fig. 1.1 Scheme for changing the propeller thrust when reversing: 
a - by formula (1.5); b - according to the thrust meter 


Thus, based on the analysis of publications on the research topic, an attempt 
was made to further improve the methodology and technique for calculating the 
inertial-braking characteristics of ships, including the following provisions: 


analytical solution of the differential equations of the slowed and acce- 
lerated movement of the ship, in which the thrust force of the propeller 
is taken equal to its average value for the work done, and the exponent in 
the motion resistance formula is equal to any integer or fractional posi- 
tive number; 

change in the thrust force of the propeller in the process of decelerated and 
accelerated movement of the ship is approximated by linear equations; 
average value of the force of resistance to the ship movement is determined 
by the theorem about its average integral value; 

theorems on the change in the momentum and kinetic energy are used as 
an alternative method for determining the characteristics of the slowed 
and accelerated movement of the ship; 
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- methodology for determining the time and speed of the ship’s accelera- 
tion to astern; 

- methodology for determining and accounting for the influence of a pass- 
ing and opponent current on the stopping distance of the ship; 

- methodology for determining and taking into account the effect of the 
ship’s stopping distance when choosing a safe speed in conditions of li- 
mited visibility. 

All calculated characteristics are compared with the results of full-scale tests 


of ships and their convergence is shown. 
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Chapter 2 

THEORETICAL JUSTIFICATION 

OF THE METHODOLOGY OF 
CALCULATING THE WAY AND TIME 
OF BRAKING OF THE SHIP 


The chapter is devoted to the development of a mathematical model, which 
forms the basis of an effective method for determining inertial-braking charac- 
teristics during active and passive braking - the most important maneuvering 
characteristics of a ship. Mathematical apparatus is also described that can be 
used by navigators to calculate: the exponent and proportionality coefficient 
of resistance to the movement of the ship; average value of the propeller thrust 
during reverse; resistance of a locked and freely rotating propeller during passive 
braking; values of the coefficient of the added mass of water during accelerated 
and decelerated movement of the ship. 

Keywords: mathematical model, propeller thrust, active and passive braking, 
added water mass, inertial-braking characteristics. 


2.1 The system of equations for the movement of the ship 
in the horizontal plane 


This system looks like this: 


dV, 


(m+Ax4) dt 





+(m+Ag WV 0=-K,-P +P -A, 
dv, 
Ca Kg) OE Neg VO =R,-P,,+A,, 
do 
a) a (2.1) 


where m — ship mass; A,, - added mass when moving along the X axis; 1, - added 
mass when moving along the Y axis; V, - projection of the ship’s speed on the 
X axis; V, — projection of the ships speed on the Y axis; w — angular ship speed; 
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J - moment of inertia of the ship relative to Z:J=m-x-R; A, — moment of inertia 
of the added masses about the Z axis; R,. - longitudinal hydrodynamic force on the 
hull; R, — transverse hydrodynamic force on the hull; P, — effective thrust force 
: the propeller: P,,. - longitudinal force of water pressure on the steering wheel; 

Py — lateral steering wheel force; A, — longitudinal aerodynamic force; A, — lateral 
aerodynamic force; Mr - moment of hydrodynamic force on the hull; M p — mo- 
ment of lateral steering wheel force; M4 - moment of aerodynamic force. 

The left-hand sides of system (2.1) contain inertial forces and moments. In the 
first two equations — the corresponding projections of the inertial force and cen- 
trifugal force, and in the third equation - the inertial moment about the vertical 
axis. In the right parts there are non-inertial forces and moments, recorded in 
general form. This system can be solved by numerical methods using a computer 
in order to simulate the movement of the ship during maneuvering. The final 
solution is possible only for special cases and under certain assumptions. 

So, the first equation of the system characterizes the movement of the ship 
along the X-axis during its acceleration and braking, therefore it allows to estimate 
the inertial-braking characteristics. The solution of the second equation, which 
describes the transverse displacement, makes it possible to obtain dependencies 
for the ship’s drift on circulation and under the influence of the wind. The third 
equation, which characterizes the angular motion, is used to assess the controlla- 
bility of ships. 


2.2 Derivation of design formulas for determining 
the stopping length and braking time of the ship when 
the propeller is reversed (active braking) 


The curvilinear trajectory described by the center of gravity of the ship 
during braking is called the braking distance. The shortest distance from the start 
of braking to the ship stop or to a given speed is called stopping distance. 

If the ship in the process of reversing does not deviate from the initial course 
and aerodynamic and other forces do not act on it (stopping distance length is 
equal to the length of the braking distance), then the law of its motion will be 
described by the first equation of system (2.1), i. e. 


(2.2) 


Fs 


(1+ 8)m —wV"- 


where k — coefficient of the added mass of water; R=uV“ - resistance to the 
movement of the ship, N; | - proportionality coefficient; a - exponent, whole 
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or fractional positive number; P,,, - average propeller thrust for the work done; 








m - ship mass, kg; dV/dt - slowdown. 

Equation (2.2) can be solved by separating the variables: 
1+k)md 
Mba (2.3) 
Po+yv" 

in which let’s introduce an auxiliary variable: 
£ 
V= (“] x, from where x = y 7 (2.4) 
(iy 


Taking into account the auxiliary variable (2.4), let’s obtain: 











1 
P\a 
(1+ 8) =| dx 4 
1+k): a 
dx = a A - ar Im Za f (alee, (2.5) 
av a P,,X P., u 
1 
h =——. 
where f(x) a 
Let's integrate equation (2.5): 
2 
es eye J Ft) (2.6) 
av u 
where C - integration constant; 
= (2.7) 


Be liege 


Let’s assume that V= Vo, x=, t=0. Let’s determine the value of the integra- 


eo 
Fe 


tion constant: 
(1+ k)m/ P,, \@ 
(“=] F(x) 
and substitute it into equation (2.6). Let’s obtain the formula for calculating the 
braking time from the initial speed Vp to its intermediate value V: 
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1 
1+k)m( P. \e 
ts Cn Ze) [F(%)-F(2)], (2.8) 
For the case of a complete stop of the ship (V=0), equation (2.8) is simplified: 
1 
t= cin Fel F(x). (2.9) 
Po KH 


To reduce equation (2.2) to a variable stopping distance length S; let’s use the 
transformation: 


dv_yav 
dt dS 
Then: 
dV 
1+k)mV —=-uV’-P.. 2.10 
Let’s divide the variables: 
1 
( tk)mVdV __ as. 
P,+uv" 


Let’s introduce the auxiliary variable x again: 











(sayn( Hey [Fey 4 (148)m/( PV" 
dS = ~ an = “( =| f(x)de, (2.11) 
av at Pix Ee u 
% 
h =——_., 
where /,(x) Le 
After integrating equation (2.11): 
2, 
1+k a 
_s+c=! in Be F(x), (2.12) 
where C, - integration constant; 
> xdx 
R(x)=|7—. (2.13) 
0 
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Assuming that at V=Vo, x=29, S=0, let’s find the value of the integration 


constant: 
2 
1+k)m( P_ \e 
6 ay Aes) 


av 


Substitute this value into equation (2.12) and obtain the formula for calculating 
the stopping distance length from the initial speed Vp to its intermediate value V: 


-(ie8m (2) [ A (%)- A(x) ]. (2.14) 


ao 


In the case of a complete stop of the ship x=0, V=0, then: 


ie Cee e . J Rls) (2.15) 


P 


av 


In most transport ships, the dependence of the hull resistance on the speed 
of movement is close to parabolic, i.e. exponent a=2 and R=uV’. Then the 
functions F(x) and F\(x9) are integrated as follows: 





= dx 
F(x,)= o =arctg x), (2.16) 
: cee , 
where arctan Xp is in radians. 
*x,dx, 1 
F(x))= | 2-2 ==In(1+ 24 (2.17) 
ee lita 2 ea) 


and, after substituting the x value, let’s obtain the formulas for calculating the 
coasting length and the braking time from the initial speed Vo to the ship stop: 


2 2 
_ (t+k)m ( s oe pater) 
P, \M) 2 Ee 2u Ti 


2k e 


0 
(1+ “im (1+ ‘ mV, 
Fi -arctg E arcs, (2.19) 
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In formulas (2.18) and (2.19), let’s introduce auxiliary variables AS and At, 
connected by functional dependencies only with the ratio of the initial resistance 
to the middle thrust: 


AS = inf 14 fs) (2.20) 
P 


av 


RB. [R 
At= |— tg |—. 2.21 
P78, Cee 


Then the stopping distance and braking time are found by the formulas: 


2 
Pom ee LAY a AS (2.22) 
or 
et k)m(Vy -V?) ae 
2R, 
when V 40. 
Pea rn At (2.23) 
or : 
Z (1+ k)m(V, -V) he 
Ry 
when V 40. 


Fig. 2.1, 2.2 show the graphs of the dependences At and AS on the ratio Ro/Pay 
for several values of the exponent a. 





0 2 4 6 8 Ro/Pav 


Fig. 2.1 Graph of the auxiliary variable A¢ for five values 
of the exponent «a» 


Chapter 2 Theoretical justification of the methodology of calculating 





AS 





0 2 4 6 8 Ro/Pav 


Fig. 2.2 Graph of the auxiliary variable AS for six values 
of the exponent «a» 


2.3 Determination of the mean value of the propeller thrust 


The average value of the propeller thrust is determined from two linear equa- 
tions: the equation for reducing the propeller thrust during reverse from the initial 
value of the ship speed Vo to the rate of change in the direction of its rotation 
speed and the equation from the rate of change in the direction of the propeller 
speed to the complete stop of the ship (Fig. 1.1). In the first case, the average pro- 
peller thrust is determined from the equation: 

VY, +V, 


P, = Ry- (Ry BJ ™. (2.24) 


In the second case, the average propeller thrust is determined from the equation: 


P,,, =P, = const. (2.25) 
Then, 
Fe i Fs. oe E3 (2.26) 


where Ro - resistance to movement at the initial speed, N; Vo — initial speed, m/s; 
Viex. - speed of the astern beginning, m/s; P, = K{pn;D* - force of the propel- 
ler thrust to reverse in mooring mode, N; K{ =0.4(H/D,)-—0.07 - coefficient 
of propeller thrust to astern; p - mass density of water, kg/m*; D, - propeller 
diameter, m. 

The pulling force (in tf) generated by the main engine can also be approxi- 


mately calculated using the following formula: 
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P.=0.01N,, (2.27) 


where N, - effective engine power, h. p. 
If N, is expressed in kW, then the pulling force (in KN) will be: 


P.=0.14N.. 


Pulling force of propellers to astern for diesel ships is approximately equal 
to 0.7 pulling force per ahead (for turbo rovers —- 0.5). 


P/ =0.7- <= 0.007N, (tf), (2.28) 
100 
P/ =0.7-0.14-N, =0.098-N, (KN). (2.29) 


2.4 Determination of the exponent value a 


The exponent a can be an integer or fractional positive number. As already 
mentioned, for most not very fast ships it is close to 2. But in the case when the 
exponent is a fractional number, it is convenient to calculate the drag value using 
decimal logarithms, i. e.: 


log R = logu + alogV. (2.30) 


The value of the exponent can be determined by plotting on logarithmic pa- 
per several points of intersection of resistance and speed as a percentage of their 
initial values. The obtained points are approximated by a straight line, taking into 
account that the ship passes about 80 % of the coasting length during the period 
of speed decrease from 100 to 50 %. 

Then the power-law dependence can be represented as a straight line accord- 
ing to the formula (2.30), in which the exponent is the slope and can be found by 
the formula: 


a= log R, log Ry (2.31) 
log V, — log V, 
Example 2.1 Fig. 2.1 shows the dependences of resistance and speed for 
a ship with a displacement of 22100 tons. Initial speeds: V,=12 knots and 
V,, =18.3 knots. 
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In Fig. 2.3, let’s remove the coordinates of the extreme points of the ap- 
proximating straight lines (dotted lines). For a speed of 12 knots, let’s obtain: 
A, (21.5; 10), B (100; 100); for a speed of 18.3 knots: A (32; 10), B (100; 100). 

Let's substitute the coordinates of these points into the formula (2.28) and 
obtain: 


— logi00-logi0 2-1 














=1.:5: 

“1 19g100—log21.5 2-133 

: _ logi00-—logi0 2-1 90 
2 Jog100-log32. 2-15 
: B 
90 ro ro toto 
so E44 HHH 
70 EEE HEE 
so 
40 
30 
20 

YT ty 

a = a 
10 0 41 394240 50 60 70 80 90 


Fig. 2.3 Dependence of resistance on speed of 12 and 18.3 knots for a ship 
with a displacement of 22100 tons 


This example also confirms the fact that with a decrease in the initial speed, 
the exponent a also decreases and approaches unity, i. e. the dependence R= f(V) 
in the region of low speeds becomes close to linear. 


2.5 Determination of the proportionality coefficient 1 


The coefficient can be determined in the following ways, most accessible in 
practice: 

1. According to the graphs of the dependence of the hull resistance on the 
speed. If there are such graphs constructed from the results of full-scale or model 
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tests, then the problem is solved relatively simply. Several corresponding values of 
resistance and speed are removed from the graph as a percentage of their initial 
values and the exponent a is calculated in the manner described in subsection 2.4. 
Then the coefficient can be determined by the formula: 


R 
ve 





W= (2.32) 


2. According to the general formula of resistance. In accordance with the 
similarity theory [1], the water resistance to the ship movement is calculated by 
the general formula: 


R= ey, from where p= oe (2.33) 


where € = 0.0085 -—0.051Fr + 0.1246F7* — total drag coefficient obtained by ap- 
proximating the impedance curve of ships [2]; Fr= v/ Jel - Froude number; 
g=9.81 m/s’; L - ship length between perpendiculars, m; Q= D* (4.854 + 
+0.492(B/ T)) — area of wetted surface, m*; D - ship displacement, t; B/T — ratio 
of the width to the average draft of the ship. 

The calculation of the wetted surface is carried out within the expected dis- 
placement of the ship. Based on the results of the calculation, a graphical depen- 
dence Q = f(T) is built, according to which, using the load scale, the value of the 
wetted surface is determined for any displacement. 

3. By the effective power of the main engine. It is known [3] that hull resis- 
tance, ship speed, towing power and main engine power (in h. p.) of a ship are 
related to each other by the following relationships: 


EPS == EPS =-n,°N,. 


Equating the right-hand sides of these equations: 


ic 
n-n,:N,= M , from where po P ae, (2.34) 





where 7 - propulsion coefficient for a laden ship is 0.70; for a ship in ballast is 0.80; 
nN, - shaft line coefficient is 0.97 if the engine room is located in the middle of the 
ship and 0.98 if the engine room is located in the stern of the ship; N, - effective 
power of the main engine, h. p.; EPS - towing power. 
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4. According to empirical formulas. The work [4] presents numerous results 
of full-scale and model tests of laden ships. 

After processing and approximating these results, a simple empirical formula 
was obtained [5]: 


R= 55 uty, form where = 55 Alt, (2.35) 


where D,,, — laden ship’s displacement, t; L - length of the ship between perpen- 
diculars, m. 
With a Froude number: 


Fr= = < 0.25, 


Jal 


which is usually the case for most transport ships, formula (2.35) gives satisfactory 
results for a laden ship. 

To recalculate the coefficient for a smaller displacement, it is possible to 
use the formula: 


1.167 
= — : 2.36 
Ht M(B ( ) 


2.6 Determination of the added mass coefficient k 


The added mass effect occurs when a rigid hull translates in a liquid. During 
the acceleration and braking of the ship, the value of the added mass of water 
turns out to be somewhat greater than during steady motion.This is due to the 
fact that at unsteady motion the hull resistance becomes somewhat higher than at 
steady-state [6], i.e. 


R=R,+AR, 


where R, - resistance at unsteady motion with an instantaneous value of speed, 
calculated as at steady motion of the ship; AR - additional resistance caused by 
the influence of the acceleration force. 

If the value of the additional resistance is combined with the inertial force, 
then the coefficient of the added mass will be equal to: 
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dV 
A, + AR: — 
oe ee (237) 
m m 


where Am - value of the added mass of water, determined by the formula [7]: 


1.5mp 
4 


Am= -T-B, (2.38) 


where T - average draft of the ship, m; B - ship width, m; m= 3.14... 
Then the mass of the ship, taking into account the coefficient of the added 
mass of water, will be equal to: 


m+Am=(1+k)m. (2.39) 


In the domestic literature [8, 9] it is accepted to consider the coefficient of 
the added mass k to be a constant value equal to 0.1, i. e. 10 % of the ship’s mass 
displacement, which is approximately suitable for medium-tonnage ships. In prac- 
tice, for large-capacity tankers and bulk carriers, the coefficient of added weight 
can be equal to 0.05 or less, and for medium-tonnage ships - 0.12 or more. Such 
errors can significantly distort the calculated characteristics of acceleration and 
braking of ships. 


2.7 An example of calculating the length and time 
of stopping distance of a ship with active braking 


Let’s take the data for the ship from article [5]. Displacement of the ship: 
22100 t; (1+) =1.1; Vo=9.15 m/s; Vyer,=6.5 m/s; Rp=688.000 H; P,}=590700 H. 

Using formulas (2.24)-(2.26), let’s determine the value of the middle thrust 
of the propeller: 


i = R, ~ (Ry — p04 Nie. so = 
= 688000 — (68800 — 590700) FS = 95378: 


P,,, = P. = 590700. 


P,, = Py, + Py, = 590700 — 73372 = 517328. 


Wy 
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Using the formula (2.18), let’s determine the stopping distance length: 


2 
g Ut kmVe (44 | _ 
2R, 


av 





1.1-22100000-9.157 688000 
= In} 1+ 


)- 1242 m. 
2-688000 517328 


Using the formula (2.19), let’s determine the stopping distance time: 
poe vee A = 
Ry Es Fs 
_ 1.1-22100000-9.15 [688000 ov [688000 = 3185=5.3 min. 
688000 517328 517328 


During acceptance tests, the following full-scale results were obtained: 
S,=1215 mand t, = 4.9 min. The ratios of the calculated and full-scale values are 
within the experimental accuracy. 





2.8 Determination of the length and time of stopping 
distance during passive braking 


Passive braking is the slow motion of the ship after the fuel (steam) is closed to 
the main engine. Up to this point, the movement of the ship was steady and occurred 
under the action of the thrust of the propellers, which was balanced by the resistance 
force. From the moment the «stop» command is given until the fuel (steam) supply 
is stopped, the ship's speed can be considered constant and equal to the speed of 
steady motion, because this period lasts only 5-10 s. The distance traveled by the 
ship during this period is defined as the product of time and speed of steady motion. 

Then slow motion begins with the main engine dumped load. There is a rapid 
drop in the propeller rotational speed to the «free rotation» mode, when the pro- 
peller, under the action of the incident flow, operates as a hydraulic turbine. Due to 
the kinetic energy of the rotating masses of the engine, it still creates a certain posi- 
tive thrust, which gradually drops to zero. The passive braking process is described 
by the following equation, in which the propeller thrust force, P;=0: 


(1+ 8)m@ suv? =o, (2.40) 
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The time and distance of passive braking are calculated from the initial speed, 
Vo, and to the final speed, V;,, which is assumed to be 0.2Vo, or to the speed of 
loss of control, whichever comes first. The solution to equation (2.40) has the 


following form: 
Balt i, 5-H ,( 4) (2.41) 
uv) \V, uu Vv, 





The propellers in the locked state can create a braking force of more than 
50 % of the hull resistance force [9, 10]. 

The value of the resistance of the locked propeller to the movement of the 
ship can be calculated using the formula from the reference book [2], after bring- 
ing it to the SI dimensions: 


2 
2 277-0; -2 (1-) 


‘ : Vv’, (2.42) 

1+ oso 2) 

D, 
in which 
477.7-D?-@-(1- BY 
L,= a (2.43) 
1+ oa 2 
D,, 


For freely rotating propellers, according to A. Kalmakov [2], the coefficient 
can be determined from the expression: 


wu, =125- Di -@(1-¥), (2.44) 


where D, - propeller diameter, m; © - disc ratio of the propeller; ‘¥ - associated 
flow coefficient; H/D, — propeller pitch ratio; V — ship speed, m/s. 
The total resistance of the hull and propeller will be equal to: 


R=pV*+u,V? at a=2, R=(u+u,)V’. (2.45) 


For approximate calculations of the resistance of a freely rotating propeller, it 
can be assumed that R,=0.2R, w+p, =1.2u at a= 2. 
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2.9 Example of calculating the length and time 
of stopping distance for passive braking 


To the data of the ship given in the example of subsection (2.6), let’s add ad- 
ditional values: propeller diameter, D,,=6.3 m; propeller pitch ratio, H/D,=1.0; 
propeller disc ratio, © =0.6; associated flow coefficient, ‘¥ =0.29; final speed 
Vp=3 m/s. 

Using the formula (2.43), let’s determine the proportionality coefficient for 
the locked propeller: 


_A77.7D20(1-¥) _ 477.7-6.3°-0.6-(1- 0.29) _ 


ny V1+0.52-17 
1+ oof ] 


4651. 








Pp 


Using the formula (2.44), let’s determine the coefficient for a freely rotating 
propeller: 


[, =125-D,?-@(1-¥) =125-6.37-0.6-0.71= 2113. 


Using formulas (2.41), lets determine the length and time of passive brak- 
ing of the ship for a locked and freely rotating propeller from the initial speed 
Vo=9.15 m/s to the final speed equal to Vz=3.0 m/s. 

For a locked propeller: 


__(1+k)m (Vv, _,)_ (t+0.1)22100000 (“48 12423527 min 
(uw, +H) (Y, (4651+8215)-9.15\ 3 





e A : (28 


)= 2007 m=11.4 cbl. 
L+H V, 4651+ 8215 3 


For a freely rotating propeller: 





_ (1+k)m [f )- eee 


] = 527s=8.8 min. 
(u, +H) LY, (21134 8215)9.15\ 3 





fe Or ) _ (1+0.1)22100000 


in( 222) = 2624 m=14.2 cbl. 
uth \Y (2113+ 8215) 


3 
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From a comparison of the results obtained, it can be seen that for a freely ro- 
tating propeller, the stopping distance of the ship is 2.8 cbl longer, and the braking 
time is 1.8 minutes longer. 

In this case, when the resistance to motion is proportional to the square 
of the ship’s speed, the path length and braking time approach their values 
asymptotically. 

Thus, as a result of the research done, calculation formulas were obtained to 
determine the characteristics of active and passive braking of the ship. Methods for 
determining the exponent and proportionality coefficient in the formula of resis- 
tance to the movement of the ship have been developed average value of the propel- 
ler thrust during reverse; resistance of a locked and freely rotating propeller during 
passive braking; values of the coefficient of the added mass of water during acce- 
lerated and decelerated movement of the ship. Examples of calculation of the length 
time of stopping distance of the ship with active and passive braking are given. 
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Chapter 3 

DETERMINATION OF THE 
CHARACTERISTICS OF THE ACTIVE 
SHIP BRAKING. ALTERNATIVE 
APPROACH 


This chapter describes an alternative approach for determining the inertia-brak- 
ing characteristics of a ship. The method under consideration compares favorably 
with those widely used in that it is convenient for use directly on the navigating 
bridge when maneuvering a ship, since it does not require experimental data - 
for calculations, it is sufficient to have general characteristics of a ship at any dis- 
placement. Also, this approach makes it possible to exclude the traditional stages 
of braking in the calculations (signal passage, passive braking and active braking) 
and to carry out the calculation continuously from the moment the signal is sent 
to the engine room to the specified speed or until the ship stops completely. 
Keywords: alternative approach, calculation method, ship maneuvering, ship 
hull resistance, braking length. 


Captains and their assistants are well aware that proper consideration of the 
inertial and braking characteristics of the ship when maneuvering in cramped 
conditions and in conditions of limited visibility is a guarantee of trouble- 
free sailing. 

The graphs of inertial-braking characteristics in the form of rectangular 
columns, recommended by IMO for practical use, do not fully reflect the charac- 
teristics of the accelerated and decelerated movement of the ship. The linear con- 
struction of these graphs makes it difficult to interpolate the results depending on 
the value of the ship’s initial speed. When maneuvering a ship at night, it is gene- 
rally impossible to use them in the darkened wheelhouse of the navigating bridge. 

In the scientific and educational literature, it is proposed to separately calcu- 
late three stages of active reversal: the time of command passage from the bridge 
to the engine room; passive braking time and active braking time. Such stages 
significantly complicate the calculation and formalization of the process; there- 
fore, such a calculation is not very suitable for practical purposes. 
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The method described below for the continuous calculation of inertial-braking 
characteristics, without dividing the process into the mentioned stages, was develo- 
ped in order to simplify the acquisition of data so important for the navigator. 


3.1 Derivation of calculation formulas 


If the ship in the process of braking does not deviate from the initial 
course (the length of the coast is equal to the length of the braking length), then 
the law of its motion is described by the differential equation: 


(+h =-y F, (3.1) 


in which the sum of the braking forces ’ F consists of the hull resistance R and 
the propeller thrust to astern P, in Newtons; k — coefficient of added mass, V - ship 
speed, m/s, m — ship mass, kg. 

The entire braking process is divided by speed into several elementary sec- 
tions, m,n+1 and assumes that in each section the work of the resistance and 
thrust forces of the propeller is equal to the work of their average values P_,, R,,,. 
Then rewrite the right-hand side of equation (3.1) as follows: 


(1+k)m S == (Rao, pa ) 


(3.2) 
In the theory of the ship, the resistance of the hull to the movement of the 
ship is usually approximated by a power dependence of the form: 


R=yV%, 


where LU — proportionality coefficient; a - exponent. 
To determine the average resistance at each section of speed, let’s apply the 
theorem on the mean of integral calculus, according to which: 


n+1 


V - u(Vv" _ Vet) 


R,, u ve V ~ Tay)? 
V ~ Wet ! (a + 1)(Vi a1) 





(3.3) 


where V, -initial speed before braking; V, 
braking section. 


- speed at the end of an elementary 


+4 
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For the most common case, when a= 2: 


Bg 5s + Vy Va Ve (3.4) 


n+1 n 


To determine the value of the aspect ratio, it is possible to use one of the 
methods described in chapter 2. Since the ship mass practically does not change 
during braking, equation (3.2) can be represented as follows: 


d{(t+kymV]__ 


dt & 


Gn, nt Wn, net ). 

Let’s multiply both sides of this equality by dt and take integrals from them. 
In this case, on the left, where the integration is in speed, the limits of the integral 
will be V, and V,,,, and on the right, where the integration is in time, the limits of 
the integrals will be 0 and ¢. As a result, let’s find that the change in the momentum 
of the ship in a given braking section is equal to the sum of the impulses of the 
forces acting on it, i. e. 


(+k), -U+k)mV «HR, tek, ot (3.5) 
Solving equation (3.5) for ¢, let’s obtain: 
— (A+k)m = 

‘ - Ra, n+ + Ly n+ Ms Yah oe 

where ¢ - braking time in seconds. For V,,, = 0: 

(1+k)mV, 
t= 7 3.7 
RK, n+t + on net 


To find the stopping distance length, let’s represent the acceleration in equa- 
tion (3.2) as: 
dv _avdS_,dV_1dv" 
dt dSdt dS 2ds— 





Then: 


(1+k)mdV* _ 


dS a IR ay +E, ). (3.8) 


mnt 


Theory and methods for calculating the inertial-braking characteristics of a ship 





Let’s multiply both sides of equality (3.8) by dS and introduce (1+k)m under 
the differential sign. Let’s obtain the expression of the theorem on the change in 
kinetic energy in differential form: 


2 
jhe) =-(Re,_, +Po,,,.)dS, (3.9) 


whence the stopping distance length at each section of the speed during braking 
is equal to: 


_ (1+k)m 2 772 
S=7R 4p Veo (3.10) 


Gn, n+ Bn, +t 


where S - stopping distance length in meters. For V,,, = 0, formula (3.10) is sim- 


n+1 
plified: 


eu (1+k)mV? 
AE Py 


Wn, nt 


(3.11) 


The calculation by formulas (3.6), (3.7) and (3.10), (3.11) will be the more 
accurate, the smaller the difference in speeds. 

Lets make one more assumption that the force of the propeller thrust is 
a linear function of the speed from the beginning of braking to the speed of the 
beginning of the reverse, V,.,,, after which it becomes constant and equal to the 
value in the mooring mode (Fig. 1.1). At a steady ahead speed, the hull resistance 
will be equal to the force of the propeller thrust corrected by the co-flow coefh- 
cient, i. e.: 


R,=(1-)P,, (3.12) 


where is ‘¥ - associated flow coefficient. 
With this in mind, let’s obtain expressions for determining the middle thrust 
of the propeller at the given speed sections: 


V, V_+V, 
a, where V,, =—" "#4, 3.13 
yo where V,, j (3.13) 


rev. 





Py = Ry - (Ry - PY) 


When V,, =V,,,, then P,, = P% i. e. propeller thrust to astern in mooring mode. 


The thrust force of the propeller to astern in the mooring mode is determi- 
ned by the formula known from the ship theory: 
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P=pKkin, Dy, (3.14) 


where p= 1020 kg/m* — mass density of sea water; K/ - ratio of the propel- 
ler thrust to astern in mooring mode; n, - rotational speed of the propeller to 
astern, rev/s; D, — propeller diameter, m. 

In the absence of information about the ratio of the propeller thrust to astern, 
its value can be determined by the formula [1]: 


Kj= 04 -0.07, (3.15) 


P 


Example 3.1 «Uelen» motor ship D,,, =22100t, (1+k)=1.1, V,=9.15 m/s, 
V,,, =6.5 m/s, R,=688.000 N, P;=590 700 N, p=8215, a =2. 

The reverse is made from full ahead to «full astern». 

Suction coefficient during acceleration to astern t=0.40 [2]. 


The results of calculating the braking characteristics are presented in Table 3.1. 


Table 3.1 


Calculation of braking characteristics of the «Uelen» motor ship 


0, i/S. Vals @,mi/s Pz,N Ria; N S,m  t,min 5 De 























9.15 8.0 8.57 559638 604959 205.8 0.4 205.8 0.4 
8.0 7.0 7.5 575730 462778 175.6 0.4 3814 0.8 
7.0 6.0 6.5 590700 347768 168.4 0.4 549.8 1.2 
6.5 5.0 - 590700 273148 = 242.7 0.7 792.5 1.9 
5.0 4.0 - 590700 = 167038 144.4 0.5 936.9 2.5 
4.3 3.0 - 590700 =: 101318 122.9 0.6 1060.0 3.0 
3:0 2.0 - 590700 = 52028 94.6 0.7 1154.4 3.7 
2.0 1.0 - 590700 19168 59.8 0.7 1214.2 4.4 
1.0 0 - 590700 2738 20.5 0.6 1234.6 5.0 


In full-scale tests, the following results were obtained: S,=1215 m and t ,=4.9 mi- 
nutes. Deviations of the calculated values from the full-scale values are within the 
limits of accuracy determined in PIIIC-89 (10 % RMS). 

Recall that the calculation of the length and stopping time with active braking of 
this ship according to the method described in Chapter 2 gave the following results: 


S=1242 m, t=5.3 min. 
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From a comparison of the results obtained, it can be seen that both methods 
are, practically, equally accurate. 


3.2 Influence of passing and opponent current 
on the braking length of the ship 


The stopping distance length during braking can be significantly influenced 
by the fair and the adverse current [1]. 

The stopping distance length on the current relative to the ground S,, will be 
calculated in stages and then summed up algebraically. 

With a following current at the first stage, let’s calculate the stopping distance 
length S' and the braking time ¢, from the initial speed to the complete stop of the 
ship relative to the water, as shown in Table 3.1. At the second stage, let’s calculate 
the path and time of the ship’s acceleration to astern to the current speed equal 
to V,, or until the ship stops relative to the ground. The equation of motion in this 
case will take the following form: 


dV, 





(14 Ay = (1-t)P/— Ry, (3.16) 
where t=0.40 — suction coefficient. 
From where 
__ G+), (3.17) 
G=)2 =k, 
and 
2 
S= ae Cam (3.18) 
2[(1-1)P’—- R,, | 
where 
2 
poh, (3.19) 





‘av 3 


The calculation is carried out according to formulas (3.17)-(3.19). 

At the third stage, the ship, under the influence of a fair current, will travel 
a path relative to the ground, equal to the product of the current speed (with 
a positive sign) by the total time of the two previous stages, i. e. 





LL = 
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(1+k)mV; 


S,.=S=S,+V,-(,+t,)=5 —{-—.—— 
gr 2 7 1 >) 2[(1—1)P’+ R,,| 


+V,+Vito, (3.20) 


exclude the braking time from formula (3.20), for which express the products 
V,t, and V,t, through the corresponding paths. Using formulas (3.17), (3.18), 
let's compose the ratio: 


nm | 


(3.21) 


St [on 


From where V,¢, = 25). 
Let’s multiply and divide by V, the product V,t,. Similarly to the previous one: 


Tg 99 
Vo Vy 


(3.22) 

Let's substitute expressions (3.21), (3.22) into formula (3.20). After a simple 
simplification, let’s obtain a formula for determining the stopping distance length 
with a passing stream: 





(3.23) 


2 
fe s(t at (1+k)mV; 


Y 2[-1)P’- R,.| 


In the case of an counter current, at the first stage, let’s also calculate the 
stopping distance length S and the braking time ¢, from the initial speed relative 
to the water V, to the stop of the ship relative to the water. At the second stage, 
let's calculate the stopping distance S, and braking time ¢, from the current 
speed V, to the stop of the ship relative to the water. The difference S—S, will 
give the stopping distance length relative to the ground. At the third stage, the ship 
under the influence of the adverse current will travel a path relative to the ground, 
equal to the current speed (with a negative sign) multiplied by the time difference 
between the two previous stages, i. e. 


(1+k)mv2 
A(R, + Py) 


‘av 


S,,=S-S,-V,(t,-t,)=S 





Vit, + Vpt. (3.24) 


At low values of the current speed, the middle thrust can be approximately 
considered equal to the propeller thrust when rotating to astern in the mooring 
mode, i.e. P,, =P’. Taking this into account, also replacing the products v,t, 
and v,t, with the corresponding stopping distances, let's obtain the formula 
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for determining the stopping distance length relative to the ground in an coun- 


ter current: 
2 
5, =S5 1-2 _ GtkymV (3.25) 
e V 2(R,,, + P’) 


0 


Example 3.2 Determine the stopping distance length for the «Uelen» motor 
ship from full speed ahead to «full speed astern» when sailing on a following 
and counter current. Current speed V,=3 knots= 1.54 m/s. 

The rest of the data are shown in Example 3.1 and Table 3.1. 

When a fair current: 


2-1.54 1.1-22100000- 2.37 








s = 1235{ 14 21793 
er 
9.15 2|0.6-590700-"219..2.37| 
When an counter current: 
5, =1235(1 ae Eanes Sean 
2{ 5907004 8219.27] 


Thus, on a fair current, the stopping distance length relative to the 
coast (ground) increased by 498 meters, and on an adverse current it decreased 
by 367 meters. 


3.3 Determining the time and speed of the ship’s 
acceleration astern 


In the practice of maneuvering, skippers often have to take into account 
the time and speed of the ship’s acceleration to astern from a stationary state to 
a stopping distance length equal to the length of one hull. 

For example, when leaving the berth without tugs, the speed of the ship to 
astern at the end of the maneuver should be high enough to minimize the effect 
of lateral drift, and at the same time safe for timely repayment of astern inertia in 
a limited water area [2]. 

Let’s make the assumption that the maximum value of the propeller thrust 
during acceleration to astern occurs simultaneously with the command and re- 
mains constant until the end of the maneuver. In this case, the differential equa- 
tion of motion will take the following form: 
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(1+8)m2 = (1-2)P -R 


a (3.26) 


‘av? 


which solve according to the method described in subsection 3.1, and obtain the 
final formulas for determining the time, t and path length, S when accelerating 
the ship to astern: 

(1+k)mV 


t=——_——__., 3.27 
(I-)P-R, de 


(1+k)mV? 


7 Tie -z] (3.28) 


where P.=Kipn;D) - propeller thrust to astern; R,,=uV?/3 - average in- 
tegral value of the hull resistance when the initial speed of the ship is zero; 
Kj= 0.4(H, /D,)-0.07 - ratio of the propeller thrust to astern; 1, - rotational 
speed of the propeller to astern, rev/s; p = 1026 — mass density of seawater, kg/m?; 
D, - propeller diameter, m; m - mass displacement of the ship, kg; k - coefficient 
of added mass, taken equal to 0.1 and (1+0.1)=1.1; V- current speed of the ship 
to astern, m/s; H — propeller pitch, m; 1. = 58(D,, /L), proportionality coefficient, 
where D,y, - ship displacement in tons. 

Let’s divide both parts of formula (3) by the length of the hull, L. By condi- 
tion, the ratio $/L is equal to one. After simple transformations, let’s obtain the 
formula for determining the speed at the end of the maneuver: 







2PL(1-7) 


ee Ta (3.29) 
(1+k)m+ ght 

Let's substitute the obtained value of speed into the formula (3.27) and obtain 
the time of acceleration of the ship to astern. 

In the case when there is no reliable information about the suction coef- 
ficient to astern, then its value can be approximately taken equal to =0.40 or 
(1-t)=0.60 during the entire acceleration time. 

Example 3.3 CSL Spirit bulk carrier. D=87585 t. L=214.97 m; m=87585000 kg; 
D,=7.240 m; H=5.185 m; H/D,=0.716; n=66 rev/m=1.1 rev/s. 


(1+k)m=1.1-87585000=96343500 kg. 


_ 58-87585 
~ 214.97 


= 23630. 
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K,=0.4-0.716-0.07=0.216. P,,=0.216-1026-1.17-7.244=738149 N. 





2-738149- 214.97 -0.6 
96343500+ = 23630 -214.97 





= 1.38 m/s = 2.7 knots. 


[ 96343500: 1.38 


93630-1387 = 310s=5.2 min. 


0.6. 738149 — 


Conclusions. Chapter 3 describes and illustrates an alternative method for 
calculating the active braking of a ship. Thanks to the use of theorems on the 
change in the momentum and kinetic energy on an elementary interval of the 
ship’s slow motion speed, it was possible to exclude the traditional stages of brak- 
ing (signal passage, passive braking and active braking) and to carry out the cal- 
culation continuously from the moment the signal was sent to the engine room to 
a given speed or to full stopping the ship. The resulting formulas made it possible 
to develop a method for determining the braking length on a following and on an 
counter current, as well as a method for determining the speed and time of the 
ship’s acceleration to astern. 
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Chapter 4 

DETERMINATION OF SHIP 
ACCELERATION AND BRAKING 
CHARACTERISTICS 


The chapter is devoted to the development of a methodology for analytical 
calculation of ship characteristics during acceleration and braking, which are 
often not represented at all in the data on the maneuvering properties of ships. 
Meanwhile, these parameters must be known for practical maneuvering to 
ensure safe navigation. For this reason, the question of finding suitable ma- 
thematical models for the movement of the ship during acceleration and brak- 
ing is very relevant. 

Keywords: analytical calculation, acceleration, braking, propeller thrust, full- 
scale test. 


Modern ideas about the maneuvering properties of ships show that the data 
on the maneuvering properties available on the ship, which it is equipped with 
in accordance with the requirements of the International Maritime Organiza- 
tion (IMO) and the normative documents of Ukraine, account for only 20 % of 
the total number of parameters that characterize the behavior of the ship while 
in motion. 

In addition, in the existing literature on this issue, dependences are given that 
do not allow obtaining the values of the time and paths traversed by ships during 
acceleration and braking with sufficient accuracy. 

Fully maneuverable properties of the ship are characterized by the data on its 
inertial qualities and controllability, which are given in [1]. In addition, large-ton- 
nage ships have maneuverability properties that differ significantly from the charac- 
teristics of medium-tonnage ships, due to the presence of large inertial forces [2]. 

Methods for calculating the characteristics of acceleration and braking are 
given mainly in the educational literature [3, 4], as well as in. 

In [4], when solving the differential equation, it was assumed that the pro- 
peller thrust during acceleration and braking remains constant, and the assigned 
number of revolutions occurs instantly, at the time of command. The last as- 
sumption can be accepted, taking into account the fact that the engine accelerates 
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in 3-5 seconds, and the process of changing the ship’s motion mode itself takes 
a considerable time. 


4.1 Theoretical substantiation and derivation 
of calculation formulas 


When the ship accelerates from the initial to the specified speed, the inertia 
force is equal to the difference between the thrust force of the propeller, corrected 
by the suction coefficient, and the force of resistance to motion, corresponding to 
the initial speed, i.e. 


(1+8)m = (1-1)P-R (4.1) 


Let’s suppose that in the initial period of time the ship is moving at a slow speed 
at a speed equal to V,, (Fig. 4.1). At point «a» the hull resistance is equal to the 
propeller thrust force, corrected by the suction coefficient. After the command «full 
speed ahead» was given, the thrust force increased to point J, resulting in an inertial 
force ab. Since in large-tonnage ships the development time of the full speed thrust 
force is an insignificant part of the time to reach the steady-state motion speed, 
in order to simplify the conclusions, let's assume that the propeller thrust force 
reaches a given value without inertia simultaneously with the command. 





Uss Uns Uf s 


Fig. 4.1 Dependence of the propeller thrust forces and resistance 
to movement on speed 


Under the influence of the inertial force, the speed of the ship will begin to 
increase. When it reaches the value of the half ahead, let’s reduce the speed of 
the propeller and balance the thrust force with the force of resistance to move- 
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ment (point d@). If to keep the propeller speed unchanged, then the balance of 
forces will come when the full speed is reached (point e). With approaching the 
speed of steady motion, the inertial force tends to zero and the ship’s speed will 
asymptotically approach its maximum value. The nature of the change in acce- 
leration from the speed of movement is shown in Fig. 4.2. At the initial moment 
of acceleration, the acceleration increases to a maximum value, then, with an in- 
crease in the speed of the translational movement, it decreases and becomes equal 
to zero with a steady movement. 

Thus, with accelerated movement, the force of the propeller thrust simul- 
taneously with the command increases to a value (1—t)P, =), and with the 
slow motion of the ship, the force of the propeller thrust simultaneously with the 
command decreases and reaches a value (1- t) P, = R, with a decrease in speed. 


dU/dt 


Us Uns Uf s 


Fig. 4.2 The nature of acceleration change 


When the ship accelerates from a given speed (or zero speed) to a steady 
speed, the value of the inertia force is equal to the difference between the thrust 
force of the propeller, P, and the hull resistance, R, corresponding to the steady 
speed, i. e. 


Grime teak (4.2) 
dt 
where  - coefficient of the added mass of the ship; dV/dt - acceleration during 
ship acceleration; t - suction coefficient, taking into account the effect of the hull 
on the propeller operation. 


To display the calculated formulas for the path and time of acceleration and 
braking, let’s express the thrust force using the formula: 


P= K,pn°D‘, (4.3) 
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in which K, - propeller thrust coefficient; p - mass density of water, kg/m’; 
n - rotational speed of the propeller, rev/s; D, - propeller diameter, m. 

The thrust ratio, K,, is a function of the relative propeller pitch A and is 
usually determined from the propeller action curves. The dependence K, = /(A) 
is close to linear and can be approximated by a linear equation (Fig. 4.3), for which 
it is necessary to compose the ratio: 


Kam = K,, —! Kinm - K, 
r A? 


o 


(4.4) 


from which we get: 
ze 
x ’ 


o 


K, = Kim —(Kinm ~~ K,,) 
where Kj,» — thrust coefficient in the mooring mode; K,, - thrust coefficient at 
steady state motion; K, - current value of the thrust coefficient; 1, A — relative 


pitch in steady motion mode and current value of relative pitch, respectively. 


Thrust coefficient 





a X Relative pitch 


Fig. 4.3 Thrust coefficient as a function of the relative propeller pitch, 
approximated by a linear relationship 


Because the: 





A V(t-@) nD, _ e nD, aeaice 
h, nD, Vy V,' 1-9 
where @ — coefficient of the associated flow, then 
Ky = Kim — [Sina Ke (4.5) 
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Substituting expressions (4.3) and (4.4) into the right-hand side of differen- 
tial equation (4.2), let’s obtain: 





(1-t)P-R= K.. [Ass Ku | -pn*Ds(1-t)-pV%, (4.6) 
where R=uV“; u - proportionality coefficient; a - exponent; V,= Vy — speed of 
steady motion. Denote: 

P= Kia pit Dit Ry =(1=t)P.=K, pr", (4.7) 
where R, - resistance of the ship's hull during steady motion; P, - propeller thrust 


at steady motion; Kj, - propeller thrust coefficient at steady motion. 
After substituting (4.7) into (4.6) and expanding the brackets: 





(1-1) P, — Ry = Pn [2% Jy uv’. (4.8) 


Let’s determine the value of the middle thrust of the propeller in a given 
range of acceleration speed: 


| a Fm — By (iv mM ) (4.9) 
V 2 


To determine the force of the middle resistance to the movement of the ship 
in the same speed range, let’s again apply the theorem on the mean of the integral 
calculus, according to which: 


- u(Vv;" = ve) 
= (a+1)(V,-V,) (4.10) 


At the most commonly used exponent, a=2, formula (4.10) takes the form: 
R,, = © (V2 +V,V,+V2). 
ao 3 1 


Let's divide the entire acceleration process into several intervals in terms of 
speed (Vj, Vo, ... V, - boundaries of the intervals) and assume that in each interval 
the work of the resistance and thrust forces of the propeller is equal to the work of 
their average values. Then rewrite the right-hand side of equation (4.2) as follows: 
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dv _ 
-R 4.11 
dt Es ‘av? ( ) 


(1+ k)— 
from which, based on the theorems about changing the momentum and kinetic 
energy, let’s obtain formulas for determining the time and length of the accelera- 
tion and braking path: 





_(i+k 

t= P. - ye Wh; (4.12) 
__(1+k)m py2 ye 
S=7p — RJ Ve |: (4.13) 


The proportionality coefficient 1 can be determined from the general resis- 
tance formula given in chapter 2: 


Remy, from where p=EPo, (4.14) 


where & = 0.0085 — 0.051 Fr + 0.1246 Fr’ - impedance coefficient of ships, obtained by 
approximating the impedance curve; Fr = v/ Vel - Froude number; g=9.81 m/s”; 
L - length of the ship between perpendiculars, m; Q= D’/°(4.854+0.492(B/T)) - 
wetted surface area, m*; D - ship displacement, t; B/T - ratio of the ship’s width 
to the mean draft. 

The ahead propeller thrust coefficient in the mooring mode can be deter- 
mined either from the curves of the action of the propellers, or by the formula: 


Kgn = 04-5 -+ 0.018. (4.15) 
D, 


The added mass effect occurs when a solid moves in a liquid. During accele- 
ration and braking of the ship, the value of the added mass can be calculated using 
the formula given in chapter 2: 


Am = SP rp (4.16) 


where T - ship draft, m; B - ship width, m; m= 3.14... 
Then the ship mass, taking into account the coefficient of added mass, will be: 


m+Am=(1+k)m. (4.17) 
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4.2 The procedure for calculating the path length 
and acceleration and braking time for «Vasily Porik» 
motor ship 


Ship displacement, D,=29100 t; length between perpendiculars, L=174 m; 
width, B=23.0 m; average draft T=9.5 m; propeller diameter, D,=5.9 m; propeller 
pitch, H=5.04 m; propeller rotation frequency for ahead motion, n,=1.415 rev/m; 
ship speed at 1, V,=6.69 m/s, full ahead speed, V/s=8.75 m/s. 

The calculation results are presented in Table 4.1: 

1. Am= josie 9.5° -23.0 = 2493079 kg; k= ee = 0.08. 

4 29100000 


2. (1+)m = 31593080 kg. 


3. W=0. 0035-05719 = = 10234; R, = 10234-6.69? = 458053 N. 
4. R,, - uh, +V/). 
5.04 


5. King = 0.4-— + 0.018 = 0.359. 
5.9 
6. P,, = 0.359-1020-1.415?-5.9* = 890134 N. 


oe so0134 (ee) - 


6.69 


_ 31593080 


31593080 


2(Pey — Reo) 


av 


[V.-V]; S= 





Table 4.2 shows the results of field measurements of the path and accelera- 
tion time of the «Vasily Porik» motor ship from a stationary state to a speed of 
12 knots (6.2 m/s). 

According to Table 4.1 and Table 4.2, the largest discrepancies between the 
full-scale and calculated curves were 5 % in time and 10 % in the path length. This 
testifies to the accuracy of the proposed method acceptable for practice. 

The calculation results are shown in Table 4.1, it is convenient to represent 
in the form of graphs dependences v = f(5S) and v= f(t), in which the accele- 
ration and braking curves asymptotically approach the ship speed during steady 
motion (Fig. 4.4, 4.5). 
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Table 4.1 
Characteristics of acceleration and braking of the «Vasily Porik» motor ship 
V;, m/s Vx, m/s Vag» M/s Se ts US zie Dhan 
Acceleration 
0 1 0.5 18.5 37.0 18.5 37.0 0.62 
1 2 1.5 61.6 41.0 80.1 78.0 1.30 
2 3 2.5 119.0 47.6 199.0 126.0 2.10 
3 4 3.5 205.6 59.0 405.0 184.0 3.00 
4 5 4.5 363.0 80.7 768.0 265.0 4.40 
5 6 5.5 774.0 141.0 1542.0 406.0 6.70 
Braking 
8.7 8 8.35 508.5 61.0 508.5 61.0 1.00 
8 7.5 7.75 543.0 70.0 1051.0 131.0 2.20 
7.5 7 7.25 985.0 136.0 2036.0 267.0 4.40 
Table 4.2 
Acceleration from motionless state of the «Vasily Porik» motor ship 
Acceleration time, Acceleration Acceleration time, Acceleration 
min-s path, m min-s path, m 

0 0 05-20 930 

02-21 186 05-55 1116 

03-18 372 06-26 1302 

04-05 558 06-58 1488 

04-45 744 
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Fig. 4.4 Dependence of speed on the path of acceleration and braking 
of the «Vasily Porik» motor ship 
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Fig. 4.5 Dependence of the speed on the acceleration and braking time 
of the «Vasily Porik» motor ship 


So, the author has developed a method for the analytical calculation of the 
characteristics of the ship’s acceleration and braking. Due to the linear approxi- 
mation of the dependence of the propeller thrust coefficient on its relative step in 
the elementary section of the acceleration and braking speed, it becomes possible 
to obtain simple calculation formulas. Comparison of the results of the full-scale 
test of the «Vasily Porik» motor ship and the results of the calculation by the pro- 
posed method shows that the relative error in the stopping distance length does 
not exceed 10 %, and the relative error in time is 5 %. 
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Chapter 5 

PRACTICAL APPLICATION 

OF INERTIAL-BRAKING 
CHARACTERISTICS OF A SHIP 


The chapter presents a universal methodology for practical calculation of the safe 
speed of a ship, taking into account its inertial-braking characteristics, which sig- 
nificantly expands the boat master’s capabilities for efficient and safe ship control. 
‘The obtained analytical dependencies make it possible to apply the automation of 
the speed selection using a personal computer directly on the navigation bridge. 
Keywords: safe speed, ship convergence, radar station, inertial-braking charac- 
teristics, braking length, navigation in ice. 


5.1 Calculation of safe speed and minimum permissible 
distance of ship convergence 


The International Regulations for the Prevention of Collisions at Sea, 1972 
(COLREGs-72) introduce such interrelated concepts as «safe speed» and «close 
proximity of ships». 

When setting a safe speed, the rules recommend, among other factors, to 
take into account the state of visibility, the maneuverability of the ship, especially 
the distance required to bring the ship to a complete stop, as well as the character- 
istics, effectiveness and limitations of radar equipment. 

COLREGs strongly recommend the use of radar equipment to detect exces- 
sive proximity and the presence of a collision hazard. In addition to qualitative 
recommendations, COLREGs-72 do not contain any quantitative assessments of 
the safe speed and excessive proximity of ships. Each master subjectively assigns 
a safe speed and determines the minimum permissible distance of approach 
of ships. Because of this, in almost every case of collisions between ships, the main 
reasons are a significant excess of the safe speed and an underestimated distance 
of the permissible approach of ships. 

In domestic and foreign literature, a large number of publications are devo- 
ted to the study of errors in technical means, in determining the elements of the 
movement of dangerous ships and assigning the minimum permissible distance 
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of convergence to them. The most complete source on this issue is the work [1, 2]. 
However, the publications do not contain the formulation and solution of the 
problem of the mutual dependence of the safe speed and the minimum permis- 
sible distance of approach of ships. 

The aim of this research is to develop a universal method for calculating the 
safe speed and the minimum allowable distance, taking into account the inertial 
and braking characteristics of the ship when using radar information [3]. 

The minimum permissible distance will be determined by the radius of the 
danger zone of ship convergence, R,,. 

When using a radar station, the radius of the danger zone of ship conver- 
gence can be determined from the following expression [4]: 


Bas. =a+25,, (5.1) 


whence the limiting value of the safe speed should be considered the speed at 
which the braking length will be equal to: 


S,=—(R,, -4), (52) 


where R,, - radius of the danger zone; 5, - length of the braking length of the 
ship at the initial speed v; a - radar parameter. 

As can be seen from formulas (5.1) and (5.2), the radius of the danger zone 
and the value of the safe speed depend on the braking characteristics of the ship 
and the radar parameter, which must take into account the technical capabilities 
of the radar to detect the danger of the echo of another ship entering the danger 
zone. The parameter a is quantitatively expressed by distance, so its value can be 
represented as a dependence: 


a=V,T, (5.3) 


where V, - relative speed of approach of ship; T - observation time required to 
detect a dangerous ship convergence with the required accuracy. 
The shortest distance of ship convergence is determined by the formula: 


a2,2a e (5.4) 


vy 


where d,, - shortest distance of ship convergence; d - detection range of a ship 
convergence; 0’ - rate of bearing change. 
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The error in determining the shortest distance of ship convergence will arise 
mainly due to inaccuracies in determining the rate of change of bearing, since 
distances are measured with a radar relatively accurately. If, in the presence of 
a danger of collision of ships, the condition d’ > d?, is satisfied, then the rela- 
tive speed of ship convergence becomes close to radial and can be measured on 
the radar screen with high accuracy. Given these circumstances, the error in the 
shortest distance will be: 


56’ 








a 5.5 
ae (5.5) 
from where: 
Vad . (5.6) 
84 


sh 


For radars operating in the centimeter wavelength range, the standard error 
in the rate of change in bearing is 60-10~° radians/s (about 0.2 deg/min). As the 
limiting error of the shortest distance, +1 mile can be taken, then the uncertainty 
of the side of the divergence of ships will be completely excluded. 

After substituting these values into formula (5.6), let’s obtain: 


“6 
Vad’ ae 0.22d?, (5.7) 
1.0 
where V, - relative speed in knots. 
The error in the rate of change of bearing is determined by the formula ob- 
tained in [4]: 


12 
8,2 = 5,2 NT” 
from where: 
T= do fe (5.8) 
dy VN 


where 6, - standard error in bearing, equal to 0.018 radians (about 1 degree); 
5,, — error in the rate of change of bearing, equal to 60 -10~° radians/s; N - number 
of echo marks during the observation time, depending on the rotation speed of 
the antenna. 
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Substituting these values into the formula (5.8), let’s obtain: 


01 
fine 00 See 1038, ea) 
60:10°¥N JN 
The rotation speed of the antennas of ship radars is not less than 15 rpm, 
therefore the value N can be obtained from the expression: 


1 
N=-T. 5.10 
i (5.10) 





Substituting expression (5.10) into formula (5.9), let’s find the minimum 
observation time for the given conditions: 


T = {4-1038? =163 s = 0.045 h. (5.11) 


Let’s express the parameter in terms of the echo signal detection range: 
a=V,-T =0.22d’-0.045=0.01d?. (5.12) 


Taking this expression into account, formula (5.1) takes the following form: 


R,, =0.01d* + 28, (5.13) 
and correspondingly, 


f 
S,= gla. -0.01d"). (5.14) 


The influence of parameter «a» on the radius of the danger zone is schema- 
tically shown in Fig. 5.1. Due to the presence of errors in the assessment of the 
collision hazard, the position of the target on the radar screen can be conven- 
tionally represented as a shaded area of a circle with a radius numerically equal 
to the parameter. Since the parameter is a function of the distance at which the 
collision hazard is determined, the echo area will decrease with decreasing this 
distance. To ensure that the echo of another ship travels from the observer ship at 
a distance of at least two stopping distances, the value of the parameter should be 
added to this distance. 

Let’s consider examples of calculating the radius of the danger zone and the 
braking length at a given speed for a specific ship, the calculation of the braking 
characteristics of which is given in chapter 3. 
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Fig. 5.1 Influence of the radar parameter on the radius of the danger zone 


Let’s represent the dependencies § = f(V) and V = f(S) (Table 3.1), when 
braking the ship to full astern, in the form of polynomials: 


S=0.07931+ 0.6263 V-0.01496V?, (5.15) 
V=0.42640.915+0.26287, (5.16) 


where S - cable; V- knots. 

Example 5.1 The ship is sailing in the open sea, away from the coast, in 
conditions of reduced visibility. Displacement of the ship is 22100 tons, speed is 
14 knots (full maneuverable ahead motion). 

Using the formula (5.15), let’s calculate the braking length when the propeller 
astern to full astern: 


S=0.07931+ 0.6263-14—0.01496-196 =5.9 cbl.=0.59 miles. 


The radar detection range of ships is 7 miles. Let’s calculate the value of the 
radius of the danger zone: 


Raz,.=0.01-72+2-0.59 = 1.67 miles. 


Example 5.2 The ship is navigating in narrowness with reduced visibility. The 
range of convergence with ships is limited to 0.9 miles, and the target detection 
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range is limited to a five-mile radar scale. Based on these conditions, it is required 
to determine a safe speed. 
Using formula (5.2), let’s find the braking length: 


S= 5(09 —0.01-5*) = 0.32 miles=3.2 cbl. 
At S,=3.2 kbt, the safe speed will be: 
V=0.426+0.9101-3.2+0.262-10.24=6.0 knots. 


If the observer’s ship, for any reason, is unable to prevent excessive approach 
to another ship or if it hears an audible fog signal ahead of its traverse, then ac- 
cording to rule 19 of the COLREGs 72 it must reduce its speed to the minimum 
until the danger of collision has passed. 


5.2 Accounting for stopping distance when ships 
navigate in ice behind the icebreaker 


To substantiate the choice of the parameters of the movement of the con- 
voy of ships in the ice behind the icebreaker, let’s designate the length of the 
convoy (from the stern of the icebreaker to the stern of the last ship) through /,, 
assuming that there are n ships in the convoy (Fig. 5.2). 


l,=d,+L,+d,+L,+..+d,+L,, (5.17) 


where d; — distance from the bow of the i-th ship to the stern of the ship in 
front (or from the stern of the icebreaker, if i=1): 


Cebreater] ——] 


yp 
tome Hy | de ge Ta a 
[= = = a 


Fig. 5.2 Scheme of the caravan formation 


The time during which the caravan travels a distance equal to its length J, 
will be equal to: 
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t, ; (5.18) 


where V; - caravan speed. 

Obviously, the sooner the caravan travels a distance equal to its length, the 
less the channel pierced by the icebreaker will drag out. 

However, it is necessary to comply with the condition that the distance 
between every two ships is not less than the stopping distance of the ship 
going behind, so that heap can be prevented in the event of a delay of the 
ship in front. 

The most favorable speed of the caravan is determined by the captain of the 
icebreaker. At the same time, it takes into account the state of the ice, meteoro- 
logical conditions and the technical condition of the ships in the caravan. Masters 
of ships in a convoy determine the braking length of their ship and maintain an 
appropriate distance between ships. 

For example, the captain of an icebreaker assigned the speed of the cara- 
van to be 6 knots. Using the formula (5.15), let’s determine the braking length 
of the ship: 


S=0.07931+ 0.6263 V—-0.01496V7= 
=0.07931+0.6263-6—0.01496-36 =3.3 cbl. 


The captain of the ship must keep at least such a distance from the icebreaker 
or the ship in front, in order, if necessary, to completely extinguish the inertia 
of the ship. 

Conclusions. As a research result, an analytical relationship is obtained 
between the safe speed and the minimum permissible distance of approach of 
ships when using radar information. The developed technique makes it possible 
to automate the choice of a safe speed and exclude the subjective factor when it is 
assigned by the navigator. 

The calculation of the braking length of the ship is shown when navigating 
in ice behind the icebreaker. 
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Chapter 6 

PRACTICAL CALCULATION 
OF INERTIAL-BRAKING 
CHARACTERISTICS 

FOR THE «ELQUI»> CONTAINER 
CARRIER 


‘The final chapter presents the result of applying the above-described method for 
determining inertial-braking characteristics using the example of the «ELQUI» 
container ship, on which the author worked as a captain. A practical calcula- 
tion, construction and implementation of working schedules of inertial-braking 
characteristics were carried out [1]. 

Keywords: laden ship, ship in ballast, active braking, passive braking, braking, 
acceleration, braking length. 


6.1 Ship data 


Ship name: «ELQUI» container ship. 
Year of construction: 1999. 

Length (full): 184.1 m. 

Length between perpendiculars: 176.0 m. 
Width: 25.3 m. 

Laden displacement: 30843 t. 

Laden draft, T,,=9.891 m. 

Ballast displacement: 17655 tons. 

Draft, T,,=4.7 m; T,=8.7 m; T,,=6.7 m. 
Main engine: Sulzer. 

Engine power: 13320 kW/17862.4 h. p. 
Propeller: right hand rotation. 

Number of blades: 5. 

Propeller diameter: 6.150 m. 

Propeller pitch: 6.101 m. 
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Stepping ratio: 0.9921. 

Disc ratio: 0.55. 

Steering wheel: semi-balanced. 
Associated flow coefficient is 0.29. 


Full speed (sea): 113 rpm. 19.5 knots (laden); 19.75 knots (ballast). 

Full speed (maneuverable): 80 rpm. 14.2 knots (laden); 14.4 knots (ballast). 
Half speed: 69 rpm. 12.1 knots (laden); 12.4 knots (ballast). 

Slow speed: 46 rpm. 8.2 knots (laden); 8.4 knots (ballast). 

Dead slow speed: 36 rpm. 6.1 knots (laden); 6.4 knots (ballast). 


Dead slow speed astern: 36 rpm. 
Slow speed astern: 46 rpm. 

Half speed astern: 69 rpm. 

Full speed astern: 80 rpm. 


Programs for calculating inertial-braking characteristics are developed in 
the Delphi programming language. The listing of programs is presented in the 
author’s thesis [2]. 


6.2 Calculation of active braking characteristics 


Initial data: 

Laden ship displacement: 30 843 t. 

Length between perpendiculars: 176.0 m. 

Width: 25.3 m. 

Propeller diameter: 6.15 m. 

Propeller pitch: 6.101 m. 

Step ratio: 0.9921 

Propeller turns to astern: 80 rpm=1.33 rpm. 

Initial speed: 14.2 knots=7.31 m/s. 

Speed before reverse: 10 knots=5.2 m/s. 

End speed: 0. 

Number of calculated points: 10. 

Ship displacement in Ballast: 17654.7 

Average draft in ballast: 6.7 m. 

Initial speed in ballast: 14.4 knots=7.42 m/s. 

The calculation results are presented in Tables 6.1, 6.2 based on them graphs 
V=/f(S) and V=/(t), which are presented below. 
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Table 6.1 
Calculation results for a laden ship 
V,, knots V,, m/s Vast, m/s Vay, m/s Sum S,cbl Sum ¢, min. 
14.2 7.31 6.58 6.94 0.64 0.28 
12.8 6.58 5.85 6.21 1.28 0.60 
11.3 5.85 5.12 5.48 1.90 0.95 
9.9 5.12 4.38 4.75 2.49 1.33 
8.5 4.38 3.65 4.02 3.02 1.74 
7.1 3.65 2.92 3.29 3.48 2.17 
Dif 2.92 2.19 2.56 3.86 2.62 
4.2 2.19 1.46 1.83 4.13 3.09 
2.8 1.46 0.73 1.09 4.30 3.58 
1.4 0.73 0 0.36 4.36 4.06 
Table 6.2 
Calculation results for a ship in ballast 
V,» knots Vi» ma/s Viet» m/s Vav, m/s Sum §,cbl Sum ¢, min. 
14.4 7.42 6.68 7.05 0.39 0.17 
12.9 6.68 5.94 6.31 0.78 0.36 
11.5 5.94 5.19 5.56 1.16 0.57 
10.1 5.19 4.52 4.82 151 0.80 
8.6 4.52 3.71 4.08 1.83 1.04 
7.2 3.71 2.97 3.34 2.11 1.29 
5.8 2.97 2.23 2.60 2.33 1.56 
4.3 2.23 1.48 1.85 2.49 1.83 
2.9 1.48 0.74 1.11 2.59 2.10 
1.4 0.74 0 0.37 2.62 2.38 





Dependence S(v) 








5, cable 


0 1 2 3 4 s 


6 e 8 
V, knots 
Fig. 6.1 Dependence of the active braking distance on the speed for a laden ship 
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Dependence T(V) 
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Fig. 6.2 Dependence of the active braking distance on the speed 
for a ship in ballast 
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Fig. 6.3 Dependence of active braking distance on speed 
for a ship in ballast 
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Fig. 6.4 Dependence of active braking time on speed 
for a ship in ballast 
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6.3 Calculation of the characteristics of the ship’s 
acceleration from a stationary state to the speed of full 
maneuvering speed 














Table 6.3 
Calculation results for a laden ship 
V,.» knots V,, m/s V2, m/s Vay, m/s Sum S,cbl Sum ¢, min. 
0 0 1.3 0.65 0.15 0.71 
2.5 1:3 2.6 1.95 0.67 1.53 
5.0 2.6 3.9 3.25 1.72 2.53 
7.6 39 5.2 4.5 3.74 3.90 
10.0 5.2 6.5 5.85 8.34 6.33 


Dependence V(S) 
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Fig. 6.5 Dependence of the acceleration path on the speed 
for the laden ship 
Dependence V(T) 
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Fig. 6.6 Dependence of acceleration time on speed 
for a laden ship 


Note. Since the dependences of the path length and the acceleration time of 


the ship asymptotically approach the value of the steady ahead speed, the calcula- 
tion should be carried out, up to approximately 0.9V,,. 
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6.4 Calculation of the characteristics of the laden ship 
braking from full maneuverability to the dead slow 
ahead speed 


























Table 6.4 
Calculation results for a laden ship 
Veenais V,, m/s V>, m/s Vaw m/s Sum S,cbl Sum ¢, min. 
14.2 731 6.92 7.12 0.67 0.29 
13.4 6.93 6.55 6.74 1.39 0.62 
12.7 6.55 6.17 6.36 2.18 1.00 
11.9 6.17 5.79 5.98 3.06 1.46 
11.2 5.79 5.40 5.60 4.04 2.00 
10.5 5.40 5.02 5.21 5.18 2.66 
9.7 5.02 4.64 4.83 6.50 3.52 
9.0 4.64 4.26 4.45 8.15 4.66 
8.3 4.26 3.88 4.07 10.38 6.35 
7.5 3.88 3.5 3.69 13.96 9.34 
Dependence V(S) 
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Fig. 6.7 Dependence of the braking distance on speed 
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Fig. 6.8 Dependences of braking time on speed 


for a laden ship 
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6.5 Calculation of the characteristics of passive braking 
from full maneuvering to stop for a laden ship and ship 
























































in ballast 
Table 6.5 
Calculation results for a laden ship 
Vi» Knots Vi» m/s Vna1, m/s Sum S, m/s Sum ¢, min. 
14.2 7.31 6.88 0.91 0.40 
13,3 6.88 6.45 1.88 0.85 
12.5 6.45 6.02 2.92 1.36 
11.7 6.02 5.59 4.03 1.95 
10.8 5.59 5.16 5.24 2.64 
10.0 5.16 4.72 6.54 3.46 
9.2 4.72 4.29 7.98 4.45 
8.3 4.29 3.86 9.56 5.65 
7.5 3.86 3.43 11.34 7.16 
6.6 3.43 3.0 13.35 9.09 
Table 6.6 
Calculation results for a ship in ballast 
V,» knots Vi» m/s Vass m/s Sum 5S, m/s Sum ¢, min. 
14.4 7,42 6.98 0.67 0.29 
13.5 6.98 6.54 1.38 0.61 
12.7 6.54 6.09 2.14 0.98 
11.8 6.09 5.65 2.95 1.41 
11.0 5.65 5.21 3.84 1.92 
10.1 5.21 4.77 4.80 2.51 
9.2 4.77 4.33 5.86 3.23 
8.4 4.33 3.88 7.03 4.11 
7.5 3.88 3.44 8.34 5/22 
6.68 3.44 3.0 9.83 6.65 
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Dependence S(v) 
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Fig. 6.9 Dependence of the stopping distance on the speed of passive braking for a laden ship 
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Fig. 6.10 Dependence of time on the speed of passive braking for a laden ship 
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Fig. 6.11 Dependence of the path length on the passive braking speed for a ship in ballast 
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Dependence T(V) 
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Fig. 6.12 Dependence of time of passive braking speed 
for a ship in ballast 


6.6 Calculation of the speed and time of acceleration 
of a laden ship at full speed astern from a stationary 
state to a stopping distance length equal to the length 
of the ship’s hull 


Initial data: 
L=176 m; D,=30843 t; Dy=6.15 m; Hy=6.101 m; Hy/Dy=0.9921; 
Ns q= 1.33 rev/s; suction coefficient =0.40; 
(1+k)m=1.1=33927300 kg; 


58 - 30843 
H=——_ 


76. 10164; K,=0.4-0.9921—0.07 =0.327; 


P,,=0.327-1026-1.337-6.154=848 983 N; 


V= a = 2.26 m/s = 4.4 knots; 
339273004 10164176 


33927300:-2.26 sae oe min, 


2 
0.6- 848983 — eee 
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6.7 Representation of inertial-braking characteristics 
in the form of analytical dependencies 


The results of calculations given in the tables of this section can be presented 
in the form of analytical dependences of the following form: S=/(V), V=/(S) and 
t=f(V), V=/(@, in which S — in cable, V- in knots, t - in minutes. These depen- 
dencies are expressed as quadratic polynomials obtained from regression analysis. 

Laden ship. Active braking from full maneuvering to full astern: 


S=—0.016806 V?+0.5574955 V—0.0877 196; 
V=0.43878735?+0.95527645+0.4329193; 
t=0.0056746V7+0.20841 19 V—0.0187392; 
V=—0.23183610?+4.3627 109¢+0.1750743. 


Ship in ballast. Active braking from full maneuvering to full astern: 


S=—0.0100063 V7+0.333229 V—0.0527337; 
V=1.23908175?+1.56807055+0.440455; 

t=0.0029449 V?+0.1251813V—0.0133596; 

V=—0.6186171t?+7.3847232t+0.1770628. 


Laden ship. Acceleration from stationary to full ahead maneuverable speed: 


S=0.1053164 V?—0.28 12967 V+0.3549369; 
V=—0.20163665?+2.84985625+0.194408; 
t=0.0416542V?+0.12444 V+0.8006972; 
V=—0.250885t?+3.5399643t—2.3622321. 


Laden ship. Braking from stationary to full ahead maneuverable speed: 


S=0.2405675 V—3.3804017 V+12.8448427; 
V=—0.03331585?+0.97635345+6.989267; 
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t=0.2215836 V-3.6018903 V+15.2008528; 
V=—0.0963039¢?+ 1.6112987¢+7.398181. 


Laden ship. Passive braking from full ahead speed to stop: 


S=0.0898451 V—0.2458343 V—-1.3230324; 
V=—0.0203645.5?+0.89195795+5.8667 285; 
t=0.1208864 V?—1.4147753 V+4.6278166; 
V=—0.0803835¢?+ 1.5927627t+6.2120366. 


Ship in ballast. Passive braking from full ahead speed to stop: 


S=0.0635215V?—0.1708933 V—-0.9439139; 
V=—0.03866175?+1.24262935+5.8803897; 
t=0.0856813 V?—1.0192293V+3.4054004; 
V=—0.1556919¢?+2.24066841+6.2405428. 


It is convenient to use the obtained analytical dependences directly on the 
navigating bridge, having previously entered the necessary polynomial coeffi- 
cients into the computer or calculator memory. 
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GENERAL CONCLUSIONS 


The performed theoretical and experimental studies makes it possible to develop 
and test in practice methods for calculating the inertial characteristics of a ship during its 
acceleration and braking, when the propeller is reversed from ahead to astern and during 
passive braking, taking into account the impact on the ship of the passing and opponent 
currents. The developed methods provide the calculation of inertial characteristics with an 
accuracy comparable to the results of full-scale tests, which makes it possible to abandon 
traditional expensive full-scale tests and at the same time significantly increase the amount 
of useful information on the navigating bridge, which increases the level of safety when 
navigating a ship in confined navigation conditions. 

‘The research results can be qualified as a theoretical generalization and practical solution 
of a major scientific problem, which is important for ensuring the safety of ship navigation. 

For the first time, a methodology for calculating the path and time of acceleration and 
braking of a ship, suitable for direct use in practice, is created and formalized. The technique is 
based on the application of theorems about the change in momentum and kinetic energy. The 
work of the forces of resistance to movement and the propeller thrust is equated, respectively, 
to the work of their mean integral and arithmetic mean values. The change in the propeller 
thrust coefficient during acceleration and braking of the ship as a function of the relative step 
is approximated by a linear relationship. Comparison of the calculated characteristics with the 
experimental ones is made for the «Vasily Porik» motor ship. The discrepancy between the ex- 
perimental and calculated data is: along the acceleration path 2 %, in the acceleration time 5 %. 

An alternative method is created for calculating and constructing working graphs of 
the path length and braking time of the ship when the propeller is reversed from ahead to 
astern. For the first time in the design formulas, the drag force to the movement of the ship 
is presented proportional to the speed value in any positive power, and the propeller stop 
force is approximated by two linear equations. The value of the exponent in the formula 
for resistance to motion is proposed to be determined as the slope of a straight line drawn 
tangent to the resistance curve plotted on a logarithmic paper. A verification calculation of 
inertial-braking characteristics for a full-scale ship is made. The calculation results are com- 
pared with the experimental data. The relative error of the calculation method is five percent. 

A method is developed for calculating the speed and time of a ship's acceleration to 
astern from a stationary state to a stopping length equal to the length of one hull. 

Software is created in the Delphi programming language. 

According to the developed methods, practical calculation, construction and im- 
plementation of working schedules of inertial-braking characteristics for the «ELQUI» 
container ship, on which the author worked as a captain, are carried out. 

Research results are versatile and effective. They can be used for transport, field, 
technical, military and other types of ships in the development of initial requirements, in 
the production of verification calculations, in the determination of inertial and braking 
characteristics of ships in service, in the analysis of accidents, in the creation of functioning 
algorithms for automatic control of processes. acceleration and braking of the ship, as well 
as in the educational process in the discipline «Ship Management». 
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